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Abstract 

This paper examines the shielding effect of the 
mean flow and refraction of sound in non-axisymmetric 
jets. A general three-dimensional ray-acoustic ap- 
proach is applied. The methodology is independent 
of the exit geometry and may account for jet spread- 
ing and transverse as well as streamwise flow gradients. 
We assume that noise is dominated by small-scale tur- 
bulence. The source correlation terms, as described 
by the acoustic analogy approach, are simplified and 
a model is proposed that relates the source strength 
to 7/2 power of turbulence kinetic energy. Local char- 
acteristics of the source such as its strength, time- or 
length-scale, convection velocity and characteristic fre- 
quency are inferred from the mean flow considerations. 
Compressible Navier Stokes equations are solved with 
a k-e turbulence model. Numerical predictions are pre- 
sented for a Mach 1.5, aspect ratio 2:1 elliptic jet. The 
predicted sound pressure level directivity demonstrates 
favorable agreement with reported data, indicating a 
relative quiet zone on the side of the major axis of the 
elliptic jet. 

Nomenclature 

a, 6 length of semi major and minor axes 

aoo freestream sound speed 

a(X) local sound speed 

C(X) normalized sound speed 

D eq area equivalent jet diameter 

k turbulent kinetic energy = vivi/2 

M Mach vector 

pi component of phase normal p 

s arc length 

T turbulence intensity =177777 /3 

U c source convection velocity 

Uj jet exit velocity 

fluctuating velocity component 
6 dissipation rate of kinetic energy 

= i '(dvi/dXjXdvi/dXj) 

#oo polar observer angle 

0 momentum thickness 

£ 1,^2 elliptic coordinates 
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p density 

r time-delay of correlation 

<f > oo azimuthal observer angle 

observer and source frequencies 

Subscripts 

o source location 

oo far field 

1. Introduction 

The aim of the present study is to determine the 
far-field directivity of sound due to a jet issued from a 
non-axisymmetric exit plane. The theoretical argument 
is a derivative of LighthilFs Acoustic Analogy, with 
modifications to account for the three-dimensional re- 
fraction effects. It is generally accepted that, in theory, 
direct numerical simulation (DNS) based on the full, 
compressible Navier-Stokes equations govern both the 
generation and propagation of sound. Current efforts 
in computational aeroacoustics 1 (CAA) are directed to- 
wards resolving various computational issues, such as 
grid resolution (for frequencies of interest), computa- 
tional domain, boundary conditions, finite difference 
scheme and so on. The presence of shock expansion 
fans, under imperfectly expanded conditions, will re- 
quire additional attention in order to predict the shock- 
associated noise. 

Apart from DNS, other computational models will 
require some sort of turbulence modeling to predict 
the mixing noise of a typical shock-free jet. For ex- 
ample, Mankbadi et. al 2,3 have opted for an extension 
of the large-eddy simulation approach (LES). This ap- 
proach, models the fine-scale turbulence and calculates 
the large-scale components from the filtered Navier- 
Stokes equations. The calculated sound field is thus 
associated with the large-scale turbulence. Further 
simplifications were proposed 4,5 by neglecting viscos- 
ity and nonlinear effects utilizing linearized Euler equa- 
tions (LEE). This last approach is less computer inten- 
sive and has been suggested as an alternative to (LES) 
for three-dimensional geometries. Sample predictions 
with LES and LEE are presented in the references cited 
above. The results are encouraging and succeed in cap- 
turing the near-field directivity of jet noise at Strouhal 
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numbers of interest. 

LighthilPs formulations, on the other hand, sug- 
gests a source correlation term based primarily on the 
unsteady velocity fluctuations. Complications associ- 
ated with numerical calculation of the LighthilPs source 
term have resulted in various modeling efforts. These 
efforts have produced alternatives to LighthilPs source 
term, such as those derived from a mean flow so- 
lution and with an appropriate turbulence modeling. 
Meanwhile, the flow-acoustic interaction is explicitly 
accounted for as an additional superimposing term. 

We start with a brief review of the source modeling 
as derived from a compressible Navier-Stokes solution 
with a k-e turbulence model (e.g., Ref. 6). Our focus, 
however, is the assessment of the mean flow effects on 
directivity of jet noise. 

Two models are proposed for acoustic predictions. 
These models differ in the application of flow-acoustic 
interaction. The first model incorporates Balsa’s 7 so- 
lution. Balsa examines the high frequency solution to 
Lilley’s equation in a cylindrical coordinate system. He 
obtains directivity factors for each of the quadrupoles 
contained within a turbulent eddy volume element. To 
extend Balsa’s solution to non-axisymmetric jets, we 
apply a quasi-3D approach. This is done by integrating 
Balsa’s solution azimuthally, i.e., treating each stream- 
wise jet slab as a representative of the entire jet volume. 
The predicted sound field, though void of any azimuthal 
directivity, may be argued to account for 3D distribu- 
tion in source strength as well as the surrounding mean 
flow. 

In the second model, a full 3D ray-acoustic ap- 
proach is applied. A system of six propagation equa- 
tions is solved numerically in a rectangular coordinate 
system. We map the CFD solution to a uniform grid, 
i.e., a grid independent of the streamwise direction, 
and apply a three-dimensional spline interpolation on 
flow parameters of interest. Mean flow gradients are 
computed readily from the ensuing spline coefficients 
as the numerical solution of refraction equations pro- 
ceeds along rays of geometric acoustics. Although a 
ray-acoustic solution succeeds in capturing the far-field 
azimuthal variation in sound, the intensive numerical 
computations suggest that the methodology should be 
considered only when the non-axisymmetric exit geom- 
etry may result in a relatively significant azimuthal di- 
rectivity. 

The source characteristic is described similarly in 
both approaches. This includes the simplification of 
fourth-order space-time correlations 6 and separation of 
space and time factors. A new function is proposed for 
the time-factor of the correlation. This function, re- 


placing the Gaussian distribution employed in the ear- 
lier models 11 , helps improve the predicted spectra. 

As an example, a shock-free Mach 1.5 elliptic jet 
with aspect ratio of 2:1 is selected. Aerodynamic pre- 
dictions are compared with data as well as predictions 
for a base convergent-divergent round nozzle, operating 
at similar conditions. The acoustic predictions include 
sound pressure level directivity (SPL), as well as 1/3- 
octave band spectra. 

2. Method of Solution 

It is generally recognized that at subsonic and 
lower supersonic speeds, small scale turbulence is the 
primary acoustic source. Each finite volume of tur- 
bulence, may be described as a multipole source that 
convects downstream and emits sound that is refracted 
by mean flow gradients. As jets becomes highly su- 
personic, large-scale structures or instability waves of 
the flow 8 become increasingly more active. In the 
present study, we assume that fine-scale turbulence is 
the source of jet noise. With a compact eddy approx- 
imation, the solution is usually expressed as a Fourier 
transform of space-time correlation function, integrated 
over the entire jet volume. Availability of unsteady flow 
solutions should, in principle, make it possible to com- 
pute correlation functions of various order and orien- 
tation. Today’s computational capabilities, perhaps, 
should encourage this effort, at least, for a selected 
number of source volume elements within the jet. The 
immediate benefit would include assessment of a num- 
ber of simplifying assumptions that have traditionally 
been employed for source modeling. For the present 
work, however, we follow an approach based on separa- 
bility of space-time correlation and rather focus on re- 
fraction and shielding of non-axisymmetric geometries. 

2.1 Source and Spectra 

As described previously, any attempt to integrate 
the source correlation function over the jet volume re- 
quires modeling of these functions. Reference 6 outlines 
a procedure used in the present study and we briefly 
discuss some of the more important steps. Roughly 
speaking, various regions of turbulence of size t , com- 
monly referred to as an eddy volume element, radi- 
ate sound to far-field that may be summed up for all 
regions of the jet if the correlation volume elements 
radiate independently. In addition, acoustical com- 
pactness which allows for the use of elementary so- 
lutions like simple sources to model complex sources 
of sound such as quadrupoles requires that ulja^ be 
small. This condition requires that the r.m.s. velocity 
fluctuation be small relative to ambient sound speed 9 
, i.e., < v 2 > 1 ^ 2 /a 00 < 1, which makes the turbulence 
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length-scale small relative to the acoustic wave length. 
It should be noted that the compactness condition is 
not that restrictive and may hold true even at moderate 
Mach numbers. Assuming that noise is dominated by 
fine-scale turbulence, contributions from the self-noise 
term, in absence of refraction and convection, is shown 
to be due to a fourth-order space-time correlation func- 
tion ViVjV r k v f £ . As is usually done, the above term is 
expressed as a linear combination of second-order cor- 
relations V{Vj. This is followed by the space-time sepa- 
ration given as u*i >'■ = Rij(()9( r )i where C and r are the 
separation vector and time-delay respectively. For ho- 
mogeneous isotropic turbulence, Batchelor 10 suggests a 
space function of the form 

Rij( 0 = T[(f + \c f)6ij - \f'CiCi/C). (1) 

Details of /(C) are provided in Ref. 10 and are not re- 
peated here. A Gaussian function is normally selected 
for the time-factor part as g(r) — exp(— r 2 /r 2 ), with 
r 0 representing the characteristic time-delay of the cor- 
relation. It was shown (Ref. 6) that r 0 may be derived 
from the ratio of kinetic energy of turbulence k and its 
dissipation rate e as \/r 0 ~ c/k. Noise spectra is calcu- 
lated from the Fourier transform of the autocorrelation 
function 


/(fi) = / g 2 (r)e ,nT dr (2) 

Jo 

which results in I(Q) = r 0 y/i r/2e ( * o) . In general 
function g(r) may be written to include other powers 
of r. For example function 

g (r) = e" c ^ T/T ^ 3 " C2 V c 3+( T / T o) 2 (3) 

with positive constants cl, c 2 and c3 transforms the 
above integral into a series of modified Bessel functions 
I\ and with added empiricism as a result of the new 
constants. With c 3 = 0 and g(r) = e -ci(T/r 0 ) 3 -c 2 |r/r 0 | 
vve have 

J(ft) = -^s/n/2Re{e r2 ercf(T)} (4a) 

2^/Ci 


where 


c 2 - M£LtqJ2) 


(46) 


and erc/(T) denotes a complementary error function. 
Constant c*i provides additional flexibility to shape the 
spectra at higher frequencies. With C 2 — 0, the solution 
reduces to that given in Ref. 11. All acoustic predic- 
tions presented here are based on equation (4a) with 
ci = 2.0 and C 2 = 1.0e-5. 


Noise spectra due to a unit volume of tu rbulence, 
and arising from the source correlation term V\Viv[v[ 
may be written as 

hm(n) = -Jj=plk?{Q.T 0 ) A ^J=Re{e r *ercf(y)} 

* (5) 
A Doppler factor relates source frequency Q to the ob- 
server frequency u>. Contributions from other source 
correlation components Ujkt are expressed similarly 6 . 

2.2 Flow- Acoustic Interaction 

It is known that the operator part of the Lilley’s 
equation accounts for refraction as well as convection. 
A close-form solution to this equation was given by 
Balsa 7 in the high frequency limit, and in a cylindrical 
coordinate system. Although Balsa’s directivity factors 
are limited to axisymmetric jets, the solution may be in- 
tegrated azimuthally for various streamwise jet slices to 
predict the SPL directivity for non-axisymmetric flows. 
This type of prediction, referred to as a quasi-3D solu- 
tion, though void of any azimuthal variation in sound, is 
very robust and should be considered a first approxima- 
tion to noise prediction for arbitrary geometries when 
far-field azimuthal directivity is relatively small. 

The general directivity of sound for jets of arbi- 
trary geometry may be assessed using ray-acoustics in 
a high frequency limit. In a high frequency approxi- 
mation, it is normally assumed that the acoustic wave 
length is shorter than the characteristic length of the 
mean flow. For high speed jets, this approximation can 
be used effectively to predict the more energetic part of 
the sound spectra and for Helmholtz numbers ( fD/a <») 
as small as one (see Tester and Morfey 12 ). When a 
parallel flow approximation is invoked, refraction of jet 
noise through a sheared mean flow of arbitrary geom- 
etry may be formulated as a two-dimensional geomet- 
ric acoustic problem 13 . It is argued 14 that for high 
speed flows, jet spreading may play a significant role in 
defining the size of cone of silence which is formed in 
the neighborhood of the downstream jet axis. In such 
cases, flow spreading may be accounted for by formu- 
lating a three-dimensional ray-acoustic problem. Thus, 
as acoustic energy, released from each finite volume of 
turbulence propagates out of the shear-layer, its spread- 
ing is inferred from the area of a ray-tube surrounding 
that ray. 

For an inviscid flow governed by linearized gas dy- 
namic equations, the coordinate X of the emitted sound 
at any point s along the ray is determined from a set 
of six equations 14,15 

=T ijPj + Mi/C i,j = 1,2,3 
as 
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dpi 

ds 



dT jk d 

dXi Pk Pl dX^ 


Mi 

c 


< 6 “> 


and 


T{j — 6{j M{Mj 


M(X) = 


U(X) 

a (X) 


C(X) = 


a (X) 


m 


Vector p is normal to the phase front and M% denotes 
the component of the Mach vector M. Equations (6) 
are solved numerically subject to initial conditions as- 
sociated with the source location and direction of emis- 
sion. Let unit vector X = ( cos /i, sin p cos 6 , sin p sin <5) 
denote the direction of emission at the source, and sub- 
scripts o and oo refer to source and observer locations 
respectively, then 


X = X 0 at s = 0 


1 M 

Po = AX + (AXM--)^, (7a) 

A = [/J 2 + (X.M)Y’/C, /? 2 = 1-|M| 2 . (76) 

As rays emerge from the shear layer, they become 
straight and the ray speed A approaches one. For every 
set of radiation angles (/j, 6) a pair of angles (#ooi <£oo) 
may be calculated from equations (6). Variation of 
the ray-tube area is then expressed as the Jacobian of 
the transformation, i.e., | d(6oo , <f> oo)/d{p, £) |- It is 
shown 14,16 that the far-field mean square pressure di- 
rectivity for a convecting quadrupole source of strength 
becomes 


P L « ( 


PoqQqo \2r>4 / QijPiPj -.2, a o 


AirR 




Ip I 2 


«5o Po 


o-S:-po) a (i-7r-pco) a f «-n/i i . 

(^^••po) 5 sine °° i a( fe , ,tr ) r 

( 8 ) 

where U c is the source convection velocity and w is the 
observer frequency related to source frequency fi as 


U = fi/(l - — • Po). 
a oo 


( 9 ) 


For an isotropic quadrupole source, QijPiPj = Q 0 | P | 2 * 
If source spectral power density Q 2 0 is replaced with 
hni/Q 4 from (5), the directivity pattern is found as 


p L*(^) 2 (p™ k ’)( QT °) 4R ^ er * erc f( r W(zt x3 opo) 

a oo 

1 


nR' 

(l-^ Po^l-^ Poo) 2 , sinp 


0-^Po) 5 


-(- 


1 0 o o | £(£°° i tz 

00 I d(i>,6) 


)■ 


(10) 


Roughly speaking, when the difference between U Q and 
U c is small, and assuming that Uco is negligible, the 
directivity pattern scales as (1 — ^ • p 0 )~ 3 outside the 
zone of silence. Near and within the boundary of zone 
of silence, the factor (^^)/| < 9 ( 0 oo> <fioo)/d(p } 6 ) | al- 
ters the above directivity and results in concentration 
of noise near the boundary of zone of silence, i.e., a 
sharp peak followed by a rapid decay into the zone of 
silence. 

2.3 Aerodynamic Predictions 

Flow predictions were made with the NPARC 
Navier-Stokes code and with a recently installed 17 low 
Reynolds number k-e model of Chien 18 . The conical 
and elliptic jets considered had an exit area of 1.571 
in 2 each, with aspect ratio of AR=2:1 for the elliptic 
jet. Both nozzles are operated at the design pressure 
ratio of 3.67 and 564° R total temperature to give exit 
Mach number of 1.5. Some aerodynamic and acoustic 
data 19 were available for validation. Because of symme- 
try, only one quarter of the elliptic jet is modeled. The 
major and minor axes are planes of symmetry modeled 
with NPARC’s slip wall boundary conditions. A grid 
having 121x91 points in axial and radial directions and 
25 points in circumferential direction was selected. The 
domain of the grid includes 10 diameters into the nozzle 
and up to 70 diameters downstream of the exit plane. 

For convenience , elliptic coordinates and £2 are 
introduced 


l/ = (icos£ 2 , z = (Zi/p)sinh, x = x, (11) 

where x denotes the streamwise direction. The above 
coordinates maintain a constant aspect ratio of p — a/b , 
with an area element dA = {£\ / p)d£\d ( ) 2 (see figure la). 
Within one quarter of the elliptic jet, the azimuthal an- 
gle 0, measured from the major axis, is equally divided 
over 25 grid points and relates to elliptic angle £2 as 
fan ^2 = ptan<j>. Two-dimensional views of the elliptic 
grid for major axis plane xy and spanwise plane yz are 
shown in figures lb and 1c. 

The centerline velocity decay for both 2D and 
3D flow predictions are sensitive to the selection of 
the spanwise boundary conditions. With a free- type 
boundary condition, the predicted centerline velocity 
decay (figure 2) shows a supersonic core length of 11.9 D 
and 11.3 D eq for round and elliptic jets compared to 
14. 2D and 11.2D C $ reported experimentally 19 . Al- 
though a slip-type boundary condition will improve the 
2D predictions (see Ref. 11), in the interest of 3D 
acoustic predictions, however, boundary conditions fa- 
vorable to elliptic jet are selected for mean flow com- 
putations. 
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Velocity profiles along the major and minor axes 
of the elliptic jet are shown in figures 3a and 3b respec- 
tively. Comparison with data of Seiner 19 demonstrates 
good agreement along the minor axis plane. The major 
axis predictions, on the other hand, appear to under- 
estimate the spreading rate. Similar conclusions were 
drawn from the axial distribution of momentum thick- 
ness defined for compressible flows 


0 = 



pu 

( pu) ct 




( 12 ) 


where subscript cl denotes the centerline values. The 
predicted momentum thickness (figure 4a) shows little 
difference along the major and minor axes near the exit. 
A more pronounced difference is observed downstream 
where thickness along the major axis shows a slightly 
faster growth rate. The round nozzle momentum thick- 
ness in figure 4a is in agreement with the elliptic jet’s 
minor axis. Data of figure 4b indicates that momen- 
tum thickness is nearly equal along both axes over the 
length of the potential core and starts to grow at a 
faster rate along the major axis near the end of the 
core. Notice that the round nozzle momentum thick- 
ness, shown here for a design Mach number of 2.0, is in 
agreement with the general behavior along the minor 
axis as was concluded from predictions of figure 4a. 


Sample prediction for Mach 1.5 elliptic jet us- 
ing a quasi-3D approach will be discussed shortly. A 
full 3D ray-acoustic solution, on the other hand, is 
obtained according to equation (10) and in conjunc- 
tion with equations (6) and (7). Here, as in the 2D 
case, source strength k 7 ^ 2 and characteristic time-scale 
r 0 are derived from CFD- ke solutions according to 
1 jr 0 — a(e/k). Empirically constant a is included 
in constant c\. Source convection velocity U c is ex- 
pressed as weighted average of exit and local velocities 
U c = 0.5U o + 0.3Uj. For an observer frequency u 
, the corresponding source frequency £1 is calculated 
from equation (9), where p c denotes the normal to the 
phase front at the source. At selected observer coordi- 
nates (i?, 0oo, <£oo), the source radiation angles (//, 6) 
and the corresponding phase normal p Q are found by 
solving (6) and (7) numerically as a boundary-value 
problem. Once the missing initial values are found, 
they may be used as a first guess for the neighboring 
source volume element AV and the process continues. 

In order to integrate the propagation equations, a 
B-spline interpolation of the CFD solutions is carried 
out in three dimensions 

N x 

Cnmi Bn t k X) t x (z) ,ti , (£l)x 

l—l m= 1 n = l 


Predicted turbulence intensity profiles along major 
and minor axes of the elliptic jet are shown in figures 
5a and 5b. Here the percent turbulence is defined as 
100 x {viVi/$) l f 2 /Uj . Both figures indicate a maximum 
level in the neighborhood of the jet lip-line and an exit 
value of 11.3%, that will gradually increase to 12.5% 
at about 5 D eq and decays farther downstream. The 
centerline value is zero within the core to about 6D eq 
and will rise to 3.4% at 7.3 D eq and 8.7% at 10.7 D eq . 
Although the existing modifications 17 in NPARC allow 
for turbulence kinetic energy input as an inflow bound- 
ary condition, the predictions within the core remain 
insensitive to the specified inflow values. Experiments 
should be performed to verify turbulence predictions. 

2.4 Acoustic Predictions 

A quasi-3D prediction model may be employed to 
predict the noise of 3D jets when the azimuthal direc- 
tivity is of less significance. This is done by integrating 
Balsa’s solution azimuthally, i.e., treating each stream- 
wise jet slab as a representative of the entire jet. The 
approach is numerically robust and should be consid- 
ered as a viable solution when azimuthal directivity is 
not a design criteria. Geometries such as round lobed- 
mixer nozzles may not result in a significant azimuthal 
directivity in the far-field. This is more so with increas- 
ing the number of lobes. 


B t ,k (3 ,t l3 (h) = f(x,Zi ,6)- (13) 

Tensor coefficients C n mi may be found for each flow pa- 
rameter of interest by solving (13) as a system of simul- 
taneous equations applied to data points (x,£i,£ 2 >/)• 
Here k X) k^ x and k^ 2 are the orders of spline , t x> t ^ 
and are the corresponding knot sequence in x, 
and £2 directions respectively 23 , and B n t k } t denotes the 
n-th B-spline of order k with respect to knot sequence 
t. It should be noted that the above process assumes 
a uniform grid, i.e., one independent of the streamwise 
direction. To this end, a postprocessing of the CFD 
solution may be necessary. This is usually done by se- 
lecting a common spanwise grid structure, say one de- 
fined at X/D eg ~ 10.0, and consequently mapping the 
flow field to this grid at each streamwise location using 
a two-dimensional B-spline interpolation. Once spline 
coefficients Cnmi in equation (13) are determined, the 
directional derivatives become readily available at any 
point within the jet. 

It is known that refraction results in a zone of si- 
lence for the high frequency noise. The zone of silence, 
denoted as 0*(X O ), and measured from the jet axis, is 
defined as the smallest polar angle that may receive 
acoustic signal from a source. This angle is found at 
each source location X Q by solving equations (6) sub- 
jet to initial conditions (//, <5) = (0, 0). Obviously, for 
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off-axis sources, a complete refraction may be achieved 
with a small value of cone angle /i if 6 is such that sound 
is directed towards the centerline. This may explain the 
observed smooth transition of SPL directivity into the 
zone of silence. However, the peak directivity level for 
a selected source is in the very neighborhood of 0* . 

Shown in figure 6 is the boundary of zone of silence 
for sources located on the major axis plane, x-y, of the 
elliptic jet. The source axial location is indicated as the 
parameter X 0 /D eq . It is clearly seen that for the most 
part, more energetic segments of the jet, i.e. sources 
near the lip-line Y 0 /a = 1.0, radiate primarily in direc- 
tion of 0* = 50°. This is consistent with the reported 
directivity for the high Strouhal number noise for Mach 
1.5 elliptic jet. Sensitivity of the zone of silence with 
respect to azimuthal source location is shown in figure 
7. Generally speaking, the flow asymmetry will result 
in an increase in the size of cone of silence as the source 
moves azimuthally closer the minor axis. 

Sound pressure level directivity for elliptic jet as 
well as Mach 1.5 conical jet are shown in figure 8. Pre- 
dictions are on a 12 foot sideline. Here, the elliptic 
jet is treated with the quasi-3D approach. This type 
of prediction, as indicated earlier, does not produce an 
azimuthal directivity. Also shown in this figure is the 
measured SPL directivity 19 . Noise measurements of 
non-axisymmetric flows, normally show lower noise on 
the deep side of a jet, such as the major axis of the 
elliptic jet. In addition, asymmetric jets usually gener- 
ate less noise compared to conical jets operating under 
similar conditions. This may be due a combination of 
factors such as turbulence level, enhanced mixing, and 
decay in supersonic core length. No data was available 
for the conical nozzle at this point, however, based on 
measurements reported at higher temperatures 24 , it is 
anticipated that the conical nozzle is at least as noisy as 
the minor axis side of the elliptic jet. Predictions of fig- 
ure 8 indicate a reasonably good agreement with data 
in the peak noise level, though the size of the zone of si- 
lence seems somewhat exaggerated. Predicted spectra 
for the two jets are shown in figures 9 and 10. 

We apply a three-dimensional geometric acoustic 
(3DGA) approach to assess the general directivity of 
sound for the elliptic jet. A matrix of five polar angles 
0^ = 50° to 110°, in increments of 15°, and four az- 
imuthal angles <j) oo = 0° to 90°, increments of 30° is 
selected. 

Shown in figures 11a and lib are the HFGA spec- 
tral predictions along the major axis side of the elliptic 
jet. Spectral measurements, converted to 1/3-octave 
band, as well as the quasi-3D predictions are shown in 
these figures. Data shows some shift in the spectral 
peak to higher frequencies with increasing polar angle 


Oqq . The prediction methods do not exhibit a similar 
shift as data. It is conjectured that this may be asso- 
ciated with simplistic replacement of the fourth-order 
space-time correlation functions. 

Azimuthal directivity (figure 12), for the most 
part, shows 2 dB lower noise along the major axis side, 
in good agreement with data of figure 8. Due to the 
immense numerical calculations associated with ray- 
acoustic predictions, the volume integration was lim- 
ited to the more energetic elements within the jet. This 
might reflect as several dB difference in the noise level 
as seen in figure 13. In addition, predicted zone of 
silence appears much larger than measurements (fig- 
ure 13). Arguably, this may relate to the absence 
of a frequency-dependent factor in the first term of 
the expansion series 15 selected for the GA predictions. 
The directivity factors suggetsted by Balsa, though ob- 
tained in a high frequency limit, contain a wave number 
factor. It is plausible that addition of the second term 
in the expansion series could improve the ray-acoustic 
predictions, in particular near the jet axis where low 
frequency noise dominates. An accurate prediction of 
the mean flow profiles and jet spreading may also have 
a significant bearing on refraction angles. For example, 
it was shown (Ref. 11 ) that an ambient Mach number 
of Moo — 0.1, will reduce the size of the zone of silence 
by more that 8 degrees. 

The significance of constant C 2 , as used in the the 
present model for the time-delay function (see Eq. 3) is 
shown in figure 14. It is seen that the predicted spectra 
show improvements relative to case when C 2 = 0.0. 

Figure 15 demonstrates the spherical directivity 
for a typical jet segment at 4.3 < X/D eq < 5.8. This 
figure suggests that the basic azimuthal pattern of jet 
noise may be captured readily by selection of appropri- 
ate source elements within acoustically active regions 
of the flow and without a need for complete volume 
integration. 

3. Concluding Remarks 

We presented a methodology for prediction of jet 
mixing noise due to noncircular jets. A shock-free el- 
liptic jet with an aspect ratio of AR=2:1 was consid- 
ered. Predictions assume that noise is generated by 
small scales of motion and is dominated by the high 
end of the spectra. 

Two prediction models were proposed. A quasi-3D 
methodology, based on Balsa’s solution, should be con- 
sidered when the azimuthal directivity is of less signifi- 
cance. This approach accounts for asymmetry in source 
distribution as well as the surrounding mean flow, how- 
ever, the predicted directivity is axisymmetric. A full 
3D ray-acoustic solution , derived from the first term in 
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the expansion series, generates an azimuthal directiv- 
ity compatible with data while the predicted zone of 
silence remains larger than measurements. 

Contributions from large-scale coherent structures 
may be accounted for via a linear inviscid stability anal- 
ysis. Stability type analysis is highly geometry depen- 
dent and has been fully explored in Ref. 22 for an 
elliptic geometry assuming a simple representation of 
the mean velocity profiles. 

As discussed, a revisit of the turbulence correlation 
function for the purpose of source modeling and valida- 
tion is appropriate at this point. Despite this needed 
improvement, ray-acoustic approach as presented here, 
sheds light on the importance of the mean flow interac- 
tion and the role it plays in shaping the directivity of 
jet noise. 
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Figure 1 b. Grid in major axis plane, 
jet exit is at X=0.0. 
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Figure 4a. Predicted momentum thickness for round 
and elliptic jets ( NPARC-ke). 
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Figure 6. Boundary of zone of silence vs. source 

location; M=1.5, M,=0.0, AR=2:1 elliptic jet. 
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Figure 7. Boundary of zone of silence vs. azimuthal 
source location. Source at £,/a=1 .0. 
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Figure 8. SPL directivity for Mach 1 .5 round 
and elliptic jets. 
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Figure 9. Predicted 1/3-Octave spectra: Mach 
1.5 round jet 
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Figure 10. Predicted 1 /3-Octave spectra; Mach 
1.5. AR=2;1 elliptic jet 
(quasi-3D predictions) 
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Figure 11 a. 1/3-Octave spectra for Mach 1.5, AR=2:1 
elliptic jet at 6_=65, <t>_=0.0. 



Figure 11b. 1/3-Octave spectra for Mach 1 .5, AR=2:1 
elliptic jet at 0_=90, 4>_=0.0. 



Figure 13. Comparison of SPL directivity; predictions 
vs. data for Mach 1 .5 AR=2:1 elliptic jet 




Figure 1 4. Predicted spectra for the elliptic jet (HFGA). 
Constant c 2 is defined in the modeling of 
the correlation function (Eq. 3). 



Figure 1 5. SPL directivity for a jet slice at 4.3<X/D <5 .8, 
within Mach 1.5, AR=2:1 elliptic jet 
(HFGA predictions) 
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